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Precision  Control  and  Randomized  Benchmarking  of  a  12-us  Class  Superconducting  Qubit 
Executive  Summary 

Superconducting  qubits  are  solid-state  artificial  atoms,  comprising  lithographically  defined  Josephson 
tunnel  junctions  and  superconducting  interconnects.  When  cooled  to  milli-Kelvin  temperatures,  these 
superconducting  circuits  exhibit  quantized  states  of  flux,  charge,  or  junction  phase  depending  on  design 
parameters.  Such  superconducting  artificial  atoms  have  already  proven  a  useful  vehicle  for  advancing 
the  scientific  community's  general  understanding  of  coherence,  quantum  mechanics,  and  atomic 
physics,  particularly  in  regimes  not  easily  accessible  with  natural  atoms  and  molecules.  Their  potential 
for  lithographic  scalability,  compatibility  with  microwave  control,  and  operability  at  nanosecond  time 
scales  make  superconducting  qubits  a  promising  candidate  for  quantum  information  science  and 
technology  applications. 

The  primary  focus  of  this  program  was  to  develop  the  calibration  techniques  required  to  eliminate  pulse 
errors.  By  eliminating  these  errors,  we  aimed  to  achieve  a  randomized  benchmarking  single-qubit 
fidelity  of  99.9%.  The  pulse  errors  we  addressed,  in  particular,  are  those  which  occur  inside  the  dilution 
refrigerator. 

To  image  the  pulse  infidelity,  we  used  the  qubit  as  an  oscilloscope  to  "visualize"  the  pulse  shape  that 
actually  arrived  at  the  qubit.  Ideally,  we  want  a  Gaussian  pulse  shape  to  arrive  at  the  qubit.  Outside  the 
dilution  refrigerator,  we  directly  measure  the  transfer  function  with  an  oscilloscope  and  use  it  then 
predistort  the  pulse  to  ensure  that  a  Gaussian  pulse  enters  the  refrigerator.  However,  inside  the 
refrigerator,  impedance  mismatch  and  frequency-dependent  loss  serves  to  distort  the  pulse.  We  used 
the  qubit  itself  to  image  the  pulse  and,  thereby,  infer  the  transfer  function  in  the  refrigerator.  We  then 
predistorted  the  pulse  to  ensure  it  arrived  at  the  qubit  as  a  Gaussian.  During  this  period  of  performance, 
we  targeted  errors  in  the  quadrature  (Q)  microwave  component  of  our  pulses.  In  a  future  program,  we 
would  like  to  extend  this  approach  to  in-phase  (I)  errors. 

Using  these  improved  pulses,  we  achieved  99.8%,  more  than  a  30%  improvement  compared  with  before 
the  start  of  this  program.  This  resulted  in  a  publication  in  PRL  [1], 

A  second  focus  of  this  program  was  to  demonstrate  improved  quantum  control  of  our  qubit  using  these 
improved  pulses,  and  this  line  of  work  led  to  three  publications  during  this  period.  In  one  work,  we  used 
the  pulses  to  implement  a  "rotary  echo,"  the  driven-evolution  version  of  the  Hahn  spin-echo  [2],  In 
another  publication,  we  showed  that  we  could  use  dynamical  decoupling  pulse  sequences  to  improve 
phase  coherence  in  coupled  two-level-systems  (in  this  case,  a  qubit  and  an  impurity)  [3],  Finally,  we  used 
microwave  pulses  to  enable  single-shot,  low-frequency  noise  measurements  in  the  milliHertz  to  100  Hz 
frequency  band  [4], 


Introduction 


Background  on  the  persistent-current  qubit 


In  this  work,  we  used  a  superconducting  persistent-current  (PC)  qubit  [5,6],  The  PC  flux  qubit  comprises 
a  superconducting  loop  interrupted  by  four  Josephson  junctions  (Fig.lA).  When  biased  with  a  static 
magnetic  flux  /dc  ~  O0/2,  where  O0  is  the  superconducting  flux  quantum,  the  system  assumes  a  double¬ 
well  potential  profile.  The  diabatic  ground  state  of  the  left  (right)  well  corresponds  to  a  persistent 
current  /q  with  clockwise  (counterclockwise)  circulation.  The  two  diabatic  energy  levels  have  an  energy 
separation  s  =  2  /q  8/dC  linear  in  the  flux  detuning  S/dC  =/dc  -  O0/2.  Higher-excited  states  of  the  double¬ 
well  potential  exist,  although  we  focus  on  the  lowest  two  states  here. 

The  two-level  system  Hamiltonian  for  the  lowest-two  states  is:  HTLS  =  —  1/2  [soz  +  A ox\.  At  detuning 
S/dc=0,  the  double-well  potential  is  symmetric  and  the  diabatic-state  energies  are  degenerate.  At  this 
"degeneracy  point,"  8  =  0  and  resonant  tunneling  opens  an  avoided  level  crossing  of  energy  A.  Here,  the 
qubit  states  are  ctx  eigenstates,  corresponding  to  symmetric  and  anti-symmetric  combinations  of 
diabatic  circulating-current  states.  At  this  point,  the  qubit  is  first-order  insensitive  to  flux  noise  (the 
energy  level  spacing  is  independent  of  flux  around  this  point),  but  maximally  sensitive  to  charge  and 
critical  current  noise.  Detuning  the  flux  away  from  this  point  tilts  the  double  well,  allowing  us  to  tune 
the  eigenstates  and  eigenenergies  of  the  artificial  atom. 

Far  from  the  degeneracy  point,  s»A,  the  qubit  states  are  approximately  ctz  eigenstates,  i.e.,  the  diabatic 
states  with  well-defined  circulating  current.  Here,  the  qubit  is  maximally  sensitive  to  flux  noise  with 
reduced  sensitivity  to  charge  and  critical  current  noise.  The  qubit  is  read  out  using  a  hysteretic  DC  SQUID 
(superconducting  quantum  interference  device),  a  sensitive  magnetometer  that  distinguishes  the  flux 
generated  by  circulating  current  states. 


Coherence  characterization  of  the  flux  qubit  (previous  AOR/LPS  program)  [7] 

In  work  funded  by  the  ARO  and  LPS,  previous  to  this  program,  we  characterized  the  coherence  of  the 
long-lived  qubit  used  in  the  present  program.  We  first  briefly  review  these  results  to  introduce  the  qubit. 
Figs.  1A-D  present  an  SEM  of  the  PC  qubit  used  in  these  measurements,  along  with  its  measured 
relaxation  (T:=12  us),  spin-echo  (T2E=23  us),  Ramsey  (T2*=2.5  us),  and  Rabi  (TRabi=13  us)  decay  traces 
when  biased  at  the  flux-insensitive  bias  point  (degeneracy  point).  At  this  bias,  the  spin  echo  is  efficient 
(T2E  ~  2T),  although  the  inhomogenous  Ramsey  is  not  (T2*  <  Ti).  Away  from  this  bias  point,  the  qubit 
becomes  more  sensitive  to  flux  noise. 

To  understand  why  the  Ramsey  time  T2*  was  less  than  the  spin-echo  time  T2E,  we  considered  the  filter 
function  picture  of  the  Carr-Purcell-Meiboom-Gill  (CPMG)  pulse  sequence  (Fig.  IE).  The  filter  function  is 
essentially  the  Fourier  transform  of  the  free-evolution  period  between  pi-pulses.  For  a  Ramsey  sequence 
(N=0),  the  free-evolution  period  is  a  "box  car"  rectangle,  which  transforms  to  a  sine  function  in 
frequency  (Fig.  IF,  N=0).  The  spin-echo  comprises  two  box  cars,  with  one  period  acquiring  a  negative 
sign  due  to  the  single  pi-pulse  which  flips  the  Bloch  vector  180  degrees  about  the  driving  field.  Its 
Fourier  transform  is  a  sinc-like  function  offset  from  zero  frequency  (Fig.  IF,  N=l).  As  more  pi-pulses  are 
added,  the  filter  function  peaks  at  higher  and  higher  frequencies.  For  noise  power  spectra  that  decrease 
with  frequency,  this  is  a  good  situation:  more  pi-pulses  lead  to  lower  integrated  noise  power  and 
therefore  longer  coherence  times. 


Fig.  1:  Coherence  characterization  of  a  long-lived  persistent-current  flux  qubit  used  in  this  work  [7].  A)  SEM  of 

the  persistent-current  flux  qubit.  B)  Relaxation  time  Tj  =  11.7  us  and  spin-echo  time  T2E  =  23  us;  C)  Ramsey 
interferometry  yielding  an  inhomogeneous  coherence  time  T2*  =  2.5  us;  and  D)  Rabi  oscillations  with  decay  time 
T Rabi  =  13  us,  all  B-D  measured  at  the  flux-insensitive  degeneracy  point.  E)  Carr-Purcell-Meiboom-Gill  (CPMG)  pulse 
sequence  contains  N  pi-pulses  that  are  aligned  with  (CP)  or  in  quadrature  to  (CPMG)  the  pi/2  pulses.  F)  The  pulse 
sequence  in  time  corresponds  to  a  filter  function  in  the  frequency  domain  that  acts  to  color  the  noise.  Increasing 
the  pulse  number  moves  the  filter  function  to  higher  frequency  and,  for  noise  power  spectra  that  decrease  in 
magnitude  with  frequency,  leads  to  longer  coherence  times.  G)  Measurement  of  the  1/e-decay  rate  as  a  function 
of  the  qubit  flux  detuning  Ob.  First-order  flux  insensitive  point  at  Ob  =  0.  Away  from  this  point,  flux  noise  increases, 
increasing  the  decay  rate.  Increasing  the  number  of  CPMG  pulses  drives  down  the  decoherence  rate  (improves 
decoherence  time).  FI)  1/e-decay  time  as  a  function  of  the  number  of  CPMG  pulses  at  Ob  =  -0.4  m  0, where  the 
qubit  is  sensitive  to  flux  noise.  Increasing  the  number  of  pulses  improves  the  decay  time  out  to  250  pulses.  Uhrig 
dynamical  decoupling  (UDD)  is  another  sequence  compared  in  this  work. 


Fig.  1G  shows  the  decay  rate  as  a  function  of  flux  (quantization  axis)  and  number  N  of  pi-pulses.  At  the 
degeneracy  point  (Ob  =  0),  the  Ramsey  decay  rate  corresponds  to  1  /  2.5  us,  while  the  decay  rate  for 
spin-echo  (N=l)  and  all  CPMG  pulses  (N>1)  is  efficient,  that  is,  1  /  2Ti.  Note  that  although  flux  noise 
dominates,  residual  charge  and  critical  current  noise  are  the  reason  the  Ramsey  decay  rate  1/  T2*  is 
greater  than  the  relaxation  rate  1/  Ti.  It  is  inconsistent  with  second-order  flux  noise  in  this  device.  Away 
from  Ob  =  0,  the  flux  noise  sensitivity  increases,  and  the  decay  rates  generally  increase.  Increasing  the 
number  of  pi-pulses  reduces  the  decay  rate  (increases  coherence  times).  In  Figure  1G,  the  coherence 
time  is  shown  to  improve  with  pulse  number  out  to  several  hundred  pulses,  in  agreement  with 
simulation. 


Summary  of  work  performed  during  the  present  funding  period 


Pulse  calibration  and  randomized  benchmarking  single-qubit  fidelity  F  =  99.8%  [1] 


In  Fig.  2A,  the  desired  Gaussian  pulse  that  we  attempt  to  apply  to  the  device  is  distorted  by  the 
frequency-domain  transfer  function  H  =  Hext  Hmt  due  to  impedance  mismatch,  mixer  imperfections,  etc. 

If  we  know  H,  we  can  invert  it  (H1)  to  predistort  the  starting  pulse  X  such  that  we  obtain  the  desired 
pulse  at  the  output,  Y  =  H(/-/_1  X)  =  X.  The  transfer  function  Hext  corresponds  to  those  distortions  that 
occur  outside  the  dilution  refrigerator,  which  can  be  found  by  measuring  the  impulse  response  of  the 
external  system  on  an  oscilloscope.  However,  there  are  also  distortions  arising  from  imperfections  inside 
the  refrigerator  which  cannot  be  calibrated  using  an  oscilloscope.  To  address  this,  we  developed  a 
scheme  whereby  the  qubit  itself  measures  the  impulse  response.  This  approach  targets  quadrature 
errors  only  (note:  in  our  new  proposal,  we  will  target  in-phase  errors).  The  resulting  improvement  in  the 
randomized  benchmarking  decay  (Fig.  2B),  and  in  the  corresponding  error  rate  versus  pulse  width 
(Fig.2C)  indicates  that  imperfections  inside  the  refrigerator  matter,  and  that  using  the  qubit  to  measure 
their  transfer  function  leads  to  a  reduced  error  rate.  Using  this  approach,  we  obtained  a  single-qubit 
randomized  benchmarking  fidelity  99.8%  for  a  4-ns  wide  pulse. 

Dynamical  decoupling  of  coupled  coherent  systems  and  the  rotary  echo  [2,3] 

We  applied  the  advanced  pulses  to  address  coherence  mitigation  in  two  experiments.  The  first  is  the 
dynamical  protection  of  an  entangled  state  of  the  qubit  and  a  two-level  system  (Fig.  2D)  in  the  space 
spanned  by  the  SWAP  operation.  We  developed  a  pulse  sequence  that  refocuses  the  coupled-system 
with  N  refocusing  pulses  (Fig.  2E),  extending  the  coherence  time  of  the  SWAP  oscillations  (Fig.  2F)  with 
the  number  of  pulses  (Fig.  2G).  This  is  essentially  a  demonstration  of  dynamical  decoupling  with  coupled 
qubits  [2], 

The  second  example  is  "rotary  echo"  (Fig.  2H),  the  driven-evolution  counterpart  to  spin-echo  [3],  With 
spin  echo,  a  pi-pulse  quickly  rotates  the  Bloch  vector  by  180  degrees  about  the  driving  axis,  and  then  the 
refocusing  occurs  during  the  qubit's  free  evolution.  In  the  rotary  echo  case,  the  qubit  is  driven  (as  in  a 
Rabi  oscillation),  and  it  is  the  driving  field  that  is  phase  shifted  by  180  degrees.  In  this  case,  refocusing 
occurs  during  the  driven  evolution.  Fig.  21  compares  the  decay  traces  of  a  conventional  Rabi  driving 
experiment  with  a  rotary  echo  experiment  for  a  particular  Rabi  frequency  (driving  amplitude),  and  Fig.  2J 
shows  the  dependence  of  the  decay  time  on  driving  field.  In  all  cases,  the  rotary  echo  outperforms 
conventional  Rabi  flopping. 

In  performing  these  measurements,  we  identified  a  new  noise  source  that  may  arise  in  qubits  coupled  to 
harmonic  oscillators.  In  this  case,  there  are  two  driving  channels:  one  is  the  desired  direct  driving  of  the 
qubit,  and  the  other  is  a  secondary  drive  as  mediated  by  the  oscillator  (in  our  case,  the  readout  SQUID) 
with  altered  amplitude  and  phase.  Fig.  2K  shows  one  consequence  of  this  effect,  in  which  despite  a 
constant  driving  amplitude,  the  observed  Rabi  frequency  changes  with  SQUID  bias  current  (oscillator 
parameters).  That  means  that  current  noise  (fluctuating  oscillator  parameters)  causes  Rabi-frequency 
fluctuations.  It  is  this  driving-field  amplitude  noise  that  the  rotary  echo  corrects  in  our  experiment. 
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Fig.  2:  Pulse  calibration  and  examples  of  advanced  dynamical  decoupling  [13 ,14].  A)  Desired  Gaussian  pulse 
becomes  distorted  by  the  transfer  function  H  due  to  imperfections  in  pulse  generation  and  transmission. 
Predistorting  with  H1  can  correct  for  these  errors.  B)  Decay  traces  for  cases  of  no  predistortion,  predistortion 
accounting  only  for  effects  external  to  the  refrigerator,  and  predistortion  using  a  transfer  function  measured  by 
the  qubit  which  incorporates  distortions  internal  to  the  refrigerator.  C)  Randomized  benchmarking  error  rate  for 
these  cases  as  function  of  pulse  width.  F  =  99.85%  is  achieved  with  a  4-ns  pulse  width.  D)  Spectrum  of  qubit 
coupled  to  a  two-level  system  (TLS).  E)  Dynamical  decoupling  sequence  that  mitigates  low-frequency  errors  in  the 
SWAP  operation  between  the  qubit  and  the  TLS.  F)  SWAP  operation  in  the  cases  of  N=0  refocusing  pulses  and  N=6 
refocusing  pulses.  G)  SWAP  decay  time  versus  number  N  of  pulses.  FI)  Conventional  Rabi  pulse  and  its 
corresponding  rotary-echo  pulse,  which  applies  a  180-degree  phase  shift  to  the  drive.  I)  Oscillation  decay  trace  for 
the  Rabi  (red)  and  rotary  echo  (blue).  J)  Decay  times  versus  Rabi  frequency  (driving  amplitude)  for  both  Rabi  and 
rotary  echo.  Rotary  echo  generally  has  a  longer  decay  time.  K)  The  readout  SQUID,  an  oscillator,  is  a  second 
channel  for  the  driving  field  to  reach  the  qubit,  modifying  its  Rabi  frequency,  shown  as  a  function  of  the  SQUID 
bias. 


Low-frequency  noise  characterization  via  repeated  Ramsey  interferometry  [4] 

In  the  present  work,  we  used  our  advanced  pulses  to  measure  the  low-frequency  noise  spectrum  in  the 
0.1  mhlz  -  100  Hz  regime  using  a  Ramsey  interferometry  technique  [4],  Low  frequency  fluctuations  in  s 
and  A  translate  into  fluctuations  in  the  Ramsey  fringes,  which  we  measured  and  analyzed  to  extract  the 
noise  power  spectral  density  (PSD).  Amazingly,  the  1/f  nature  is  consistent  with  our  previous  work 
(previously  funded  by  ARO/LPS)  across  10  decades  in  frequency.  If  we  extend  the  low-frequency  flux 
noise  out  to  the  qubit  frequency  A  =  5.4  GFIz,  it  matches  the  strength  required  for  T:  relaxation  noise 
(the  flux  noise  couples  transversally  only  at  this  frequency).  This  suggests  that  the  low-frequency  flux 


noise  responsible  for  dephasing  may  be  the  same  mechanism  for  relaxation,  motivating  us  to  measure 
the  noise  at  higher  frequencies  (future  work). 
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Fig.  3:  Characterizing  low  frequency  noise  using  advanced  pulse  sequencing  techniques  [4].  A)  Noise  power 
spectral  density  (PSD)  for  both  ctx  (5A)  and  az  (5e)  noise.  In  the  present  work,  we  used  the  advanced  pulses  to 
measure  low-frequency  noise  (portion  of  plot  at  frequencies  f<  1 00  Hz)  using  a  repeated  Ramsey  interferometry 
technique  [4].  The  other  portions  of  this  plot  were  performed  under  a  previous  ARO/LPS  program  [7]. 
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